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8.  Cooplete 


OF  )£T20R0L00IC.'I.  FIELDS 
A.  S,  Zionin  and  A.  U.  Obukhov 

<  On  the  basis  of  the  solution  of  the  problea  of  slight  vari¬ 
ations  in  a  baroclinic  atmosphere,  a  general  classification 
of  the  main  types  of  dynamic  processes  in  the  atmosphere  ’ 

(horizontal  vorticity  motions  and  gravitational  and  acoustic  .. 
waves)  is  given.  The  work  also  gives  the  general  form  of 
the  invariant,  with  which  the  critical  stationary  state  of 
the  atmosphere  can  be  computed  using  arbitrary  initial  data, 
without  analysing  the  wave  processes  which  cause  reorgard- 
'  sation  of  the  fields.  The  "filtering"  role  of  the  quasi¬ 
static  approximation  is  explained t  it  "filters  out"  internal 
acoustic  waves  end  s oner hat  overrates  the  frequencies  of 
gravitatlonal'waves.  In  particular,  it  la  shown  that  only 
several  minutes  are  required  to  establish  quasiststic  equi¬ 
librium  in  the  atmosphere,  (  j 

■  <v  • 

As  is  known,  the  coapleto  system  of  the  equations  of  atmospheric  hydro-  ‘ 
dynamics  is  a  time  system  of  the  fifth-order  (it  contains  derivatives  of  three 
components  of  velocity,  pressure,  end  density  with  respect  to  time).  Accord¬ 
ingly,  in  order  to  solve  the  Cauchy  problem  for  this  system  of  equations  one 
must  know,  at  the  initial  moment  of  time,  the  Helds  of  five  meteorological 
elements.  With  arbitrary  initial  data,  motions  developing  in  the  atmosphere 
can  be  separated  into  relatively  slow  ones  (synoptic)  and  rapid  ones  (wave). 

The  consideration  of  rapid  wave  motions  essentially  complicates  the  analysis 
of  synoptic  processes.  Therefore,  when  investigating  synoptic  processes,  the 
equations  of  atmospheric  hydrodynamics  ere  simplified  so  that  the  simplified 
equations  will  give  an  adequate  enough  description  of  the  synoptio  processes, 
but  dll  not  contain  rapid  wave  motions  among  their  solutions.  .Such  approxi¬ 
mations  (quasistattc  and  qua  si  geos  trophic  approximations)  lead  to  dynmaio 
equations  of  the  first  timo  order,  so  that  to  describe  the  evolution  of  synop¬ 
tic  processes  one  need  know  the  Initial  field  of  only  one  meteorological  ele¬ 
ment  (e.g.,  pressure  or  a  stream  function).  The  order  of  the  dynamics  equa¬ 
tions  is  reduced  (from  fifth  to  first)  through  the  exclusion  -  "filtering  out" - 
of  certain  solutions  of  the  Initial  Bystem  of  equations  (aotually,  rapid  wav* 
motions)  hy  the  above-mentioned  simplifications.  To  prove  the  validity  of  such 
simplifications  one  must  explain  the  nature  of  rapid  wave  motions. 

Under  actual  conditions,  rapid  wave  processes  occur  with  smell  amplitudes, 
i.e.,  they  are  of  the  nature  of  slight  variations.  Accordingly,  verification 
of  tho  quasi  static  and  quasi  geos  trophic  approximations  requires  study  of  the 
nature  of  slight  atmospheric  venations;  the  following  results  ere  obtained. 
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Raoid  wav*  notions  occur  if  static  and  geos  trophic  equilibrium  ia  dis¬ 
rupted  in  a  certain  region  of  the  atmosphere.  The  resultant  waves  ere  scat¬ 
tered,  due  to  which  the  meteorological  fields  adjust  to  one  another  so  that  tha 
atsBsphore  approaches  the  state  of  static  and  geos  trophic  equilibrium.  This 
process  is  called  the  adaptation  of  meteorological  fields.  In  the  linear 
theory,  the  state  of  the  statio  end  geostrophic  equilibrium  is  stationary  and 
in  order  to  distinguish  this  stationary  state  with  arbitrary  initial  data  we 
need  only  filter  out  the  rapid  wave  motions. 

Tha  evolution  of  the  adapted  fields  (i.e.,  those  close  to  the  state  of 
static  and  geostrophic  equilibrium)  occurs  as  the  result  of  non-linear  effects, 
viz.,  absolute  vorticity  transfer  and  entropy.  These  non-linear  effect  contin¬ 
ually  create  a  tendency  toward  the  destruction  of  static  and  geostrophic  equi¬ 
librium.  But  here,  due  to  the  generation  end  scattering  of  rapid  waves, 
Captation  continually  occurs,  and  the  meteorological  fields  remain  close  to 
the  state  of  statio  aid  geostrophic  equilibrium  f adapted).  The  evolution  of 
the  adapted  fields  is  also,  strictly  speaking,  a  synontio  process  and  whan 
describing  it,  one  must  consider  its  cause  -  the  non-linear  effects  -  but  in 
place  of  a  description  of  oontinually  occurring  uiaptation  one  auy  substitute 
the  requirement  that  the  meteorological  fields  remain  adapted  at  all  times. 

This  requirement  also  provides  the  basis  for  the  quasi  static  and  quasigeo- 
s trophic  approximations. 

A.  V.  Obukhov  (1]  first  suggested  the  adaptation  of  meteorological  fields 
as  a  method  of  verifying  the  quasi  geostrophic  approximation.  In  so  doing,  ho 
sxaained  the  ease  of  a  barotropic  atnosphere.  Charaey  [2]  mode  the  preliminary 
analysis  of  the  "filtering  role"  of  the  quasistatic  and  quasigeostrophio  ap¬ 
proximations.  I.  A.  Kibel*  [3]  and  A.  S.  lonin  [It]  described  the  processes  of 
adaptation  in  a  quasistatic  baroclinic  atmosphere.  In  this  latter  work  an 
Invariant  was  found  which  is  a  generalisation  of  the  "potential  vorticity"  for 
a  quasistatic  baroclinic  atmosphere  and  it  was  established  that  in  addition  to 
the  internal  gravitational  waves  investigated  in  [3],  in  *  baroclinic  atmos¬ 
phere  more  rapid  "two-dimensional"  gravitational  waves  may  occur,  which  arc 
analo(pus  to  the  wavea  in  a  barotropic  medium. 

The  question  of  tho  adaptation  of  meteorological  fields  to  the  state  of 
quasistatic  equilibrium  and  of  the  filtering  role  of  the  quosistatlc  approxi¬ 
mation  had  not  been  studied  in  sufficient  detail  until  now.  The  present  work  - 
aims  to  fulfill  this  need  and  to  give  tha  most  general  formulntion  of  tha ' 
problem  of  tho  adaptation  of  meteorological  fields  in  a  btroclinic  atmosphere. 


The  study  of  slight  atmospheric  variations  is  significant  not  only  for 
describing  the  process  of  adaptation,  but  it  is  of  independent  Interest  In  a 
number  of  problems  of  atmospheric  physics,  e.g.,  a  certain  family  of  slight 
variations  is  the  basic  object  of  study  in  atmospheric  acoustics.  AoousUo 
variations  are  of  interest  for  atmospheric  physics  since  acoustic  waves,  which 
form  in  the  troposphere,  penetrate  into  the  upper  layers  of  the  atmosphere  end 
are  absorbed  there,  can  transmit  energy  from  the  loeer  to  the  upper  atmosphere. 
Internal  gravitational  waves  arias,  e.g.,  during  the  atreanlining  of  mount sine 
by  an  air  etrean  (aee  the  woria  of  lyre  [6}  end  A.  A.  Dorodnitsyn  [$]). 

1.  THE  EQUATIONS  0?  SLIGHT  ATI  ©SPHERIC  VARIATIONS 

ne  will  examine  the  atmosphere  as  e  liquid  in  which  processes  occur  poly* 
tropically  with  the  exponent  of  polytropy  x.  In  other  words,  we  will  consider 
that  when  liquid  particles  move,  the  magnitude  pp~V*  is  preserved  in  them, 
where  p  and  pare  pressure  and  density.  In  s  real  atmosphere  (outside  the 
bounding  layers)  the  processes  occur  adiabatically,  so  that  x  *  Cp/oy,  which  is 
equal  to  the  ratio  of  the  specific  heat  of  air  at  constant  pressure  end  constat . 
volume.  However,  when  interpreting  certain  formulas  it  will  be  convenient  if 
we  can  examine  polytropic  processes  with  « r*  arbitrary  exponent. 

As  the  basic  aquations  of  atmospheric  dynamics,  lot  us  examine  the  equa¬ 
tions  of  motion  (disregarding  dissipative  factors*)),  the  equ;  tlcn  of  continu¬ 
ity,  and  the  equation  of  the  polytropic  procoes 

•  >$-- %+¥. 

(i) 

_ _ 

«) 

There  would  be  no  great  difficulties  in'  computing  the  dissipative  forces, 
but  the  computations  would  become  more  unwieldy.  When  computing  the  dissipative 
forces  in  the  form  of  the  components  ypiv,  and  vp6r  in  the  right-hairi 

parts  of  tha  first  three  equations  of  (1)  (v  is  the  visooslty  coefficient),  a 
basic  change  in  the  results  given  below  would  result  in  the  ^pearance  of 
- 

damping  factors  a  in  tha  formulas  for  tha  amplitudes  of  harmonic  waves 
with  wave  nuaber  k. 


-  4  - 


Hor«  t  is  ttno,  u,  »,  and  w  are  the  velocity  components  along  the  axes  of 


the  cartesian  coordinates  x,  y,  and  a  (the  z-axis  is  directed  vertically  up- 

i 


ward).  —  -  —  ♦  u—  is  the  synbol  of  the  individual  derivative. 

<u  dt  dx  dy  da  * 

e  the  acceleration  of  gravity  and  t  the  Coriolis  -parar.eter.  from  now  on  to  bo 

considered  a  constant 

Tfe  nill  stutty  slight  variations  which  sight  arise  in  the  atmosphere 
against  the  background  of  one  of  its  basic  states  (whose  characteristics  will 
bo  designated  by  a  line  aver  the  letters).  The  features  of  the  slight  vari¬ 
ations  will,  naturally,  depend  on  the  properties  of  the  basic  state.  In  the 
present  work  we  will  not  explain  how  atmospheric  notion  in  its  basle  state 
affects  the  nature  of  alight  variations  and,  in  connection  with  this,  we  will 
use  as  the  basic  state  the  state  of  rest  where  pressure  p  and  density  p  depend 
only  on  s  and  are  connected  by  the  statics  equation 

g--<P  (2) 

The  relationship  between  slight  variations  and  properties  of  basic  atmos¬ 
pheric  notion  can  be  explained  in  a  separate  paper.  The  conclusions  of  the 
present  work  relative  to  emall-scalo  waves  are  also  applicable  in  presence  of 
basic  atmospheric  notions  (provided  the  velocity  field  in  the  basic  state  does 
not  contain  small-scale  lnhocogeneities). 

The  distinjjiiahing  characteristics  of  the  basic  state  of  the  ataosphere 
is  the  temperature  profile  7(a),  instead  of  which  it  will  be  convenient  to  use 
the  magnitude 

o2  -  -  niff  (3) 

P 

which  has  the  sense  of  the  square  of  the  speed  of  sound,  let  us  also  introduce 
a  soecial  designation  for  the  parameter  of  theraal  stability 


Changes  of  t  with  latitude  make  possible  tho  appearance  in  the  atrosohere 
of  additional  inertial  variations  -  the  so-called  Rossby  waves.  Those  waves 
are  oC  interest  only  when  analyzing  processes  on  a  synootic  scale;  ccmov.ting 
them  would  introduce  no  essential  changes  into  the  results  of  the  prosent  work 
on  gravitational  and  acoustie  variations  -  they  would  only  cake  the  computations 
core  unwieldy.  A  sufficiently  detailed  description  of  Rossby  waves  is  given, 
e.e.,  in  A.  U.  Xaglca's  survey  article  (7). 


(u-i)e  ♦ 


Her*  ft  Is  the  gss  constant  and  Ya  and  Y  arc  the  adiabatic  and  the  actual, 
vertical  temperature  gradients. 

Unearlsing  equations  (1)  with  respect  to  the  state  of  rest  ve  mill  get 
the  equations  of  slight  atmospheric  variations  in  the  fens 

4— •£+<*  • 

i»»  ^  _  it  _ 

(S) 

a»  /  *?■  ,  it  .  i»\ 

:  *  ““l"3r+  »  +  *  r 

where  u,  v,  w,  p  and  p  are  the  characteristics  of  slight  variations  (let  us 
assume  for  them  the  same  designations  as  for  the  initial  meteorological  ele¬ 
ments).  The  parameter  X,  equal  to  unity  in  the  general  case,  is  Introduced 
here  so  that  in  the  future  it  will  be  convenient  to  trace  its  influence  on  the 
alight  variations  of  tha  vertical  accelerations  of  particles.  Actually,  if  we 
do  not  compute  vertical  accelerations,  i.e.,  if  we  describe  the  slight'varl- 
stions  in  the  quasistatic  annroximation,  ne  should  assume  that  X-»0. 

There  are  no  difficulties  in  obtaining  the  linearised  equations  (5)  from 
tha  initial  equations  (l)  except,  perhaps,  the  last  of  these  equations.  Direct 
linearisation  of  ths  last  equation  of  (1)  yields 

Using  equations  (2)  -  (U)  as  characteristics  of  the  basic  state  of  the 
atmosphere,  let  us  transfora  the  expression  in  parentheses  in  the  right-hand 
part  of  this  last  equation  as  follows* 

Substituting  this  value  in  (6)  and  eliminating  i^t,  we  get  the  last  equa¬ 
tion  of  (S)  by  using  tha  fourth  equation  of  eystsa  (5). 

In  the  future  it  will  be  convenient  to  Introduce  the  new  unknown  magni¬ 
tudes  *,  f,  and  x  instead  of  the  velocity  components  u,  v  and.  w,  setting 


The  relationship* 


<?) 


H'-v  +  iJ*  rw-x- 
+  &-£-**•  • 


(8) 


a2  a2 

where  A  »  — x*  — k  is  the  symbol  for  the  two-dimensional  Laplace  operator. 

ax2  ay 

•how  that  the  velocity  potential  tp  and  the  stream  function  f  are  determined  by 
plane  divergence  and  the  vertical  oomponent  of  the  relative  vorticity  of  the 
velocity  field,  resoectively.  The  .magnitude  y  has  the  sense  of  the  vertical 
stream  of  mass. 

Differentiating  the  first  equation  of  system  (5)  with  respect  to  x  and 
the  second  with  respect  to  y  and  adding  the  results  no  get 

A  -  lAifr  -  Ap  . 

oXt 


Analogously,  differentiating  the  second  equation  of  system  (5)  with  re¬ 
spect  to  x  and  subtracting  the  first  aquation,  differentiated  with  respect  to 
7,  get 


-<A<? 


Under  the  natural  requirement  that  functions  <p,  and  p  be  regular  to 
Infinity,  the  sign  of  the  Laplace  operator  in  the  last  Wo  equations  can  be 
eliminated.  In  addition,  using  the  designations  of  (7)  in  the  last  three 
equations  of  (5)  we  mar  rewrite  equations  (5)  in  ths  fora 

it.  --.c.A?_fcc -.«•£,  (?) 

Ths  squstions  of  slight  vsrlstions  will  be  most  convenient  in  this  form 
for  subsequent  analysis.  v«'e  will  study  the  solution  of  these  squstions  using 
srbitrsiy  initial  data.  ’Then  . 

*"0:  *  -  V.  0;  ?- 1.(*.  if.  y,«); 

t  •  ■  M*.  if.  *);  p  -  />.(*.  y,  *), 
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Sine*  the  system  of  equation'  (9)  contains  derivatives  with  respect  to  •  • 
and  * s  of  the  sc  vid  z-order,  two  boundary  conditions  with  respect  to  s  must  bo 
formulated  to  solve  it.  The  natural  boundary  conditions  are  - 

x-»0(i-»0),  x~*0(i-oo).  (H) 

The  first  of  theso  conditions  is  obvious.  The  second  indicates  that  the 
vertical  stream  of  boss  should  revert  to  zero  at  the  lower  boundary  of  the 
atmosphere*). 

'  2.  THE  STATION/Jtt  SOLUTION  S®  THE  INVARIANT  . 


Equations  (9)  have  a  certain  family  of  stationazy  solutione  t(i  ^  p#, 
and  p#.  For  each  of  the  stationary  solutions  the  following  relationships 

ere  fulfilled  '  t  H 

f,mO,  x,mO,  p.eslfc,  J,™— y-57-.  (12) 

so  that  stationary  notions  are  horizontal  and  non-d ivergent  and  the  formulas 
of  geos  trophic  wind  and  the  statics  equations 

U.  jr  »i--t -|J“  ■  “j“ “  ~ It,’  VU/ 

are  applicable  to  then. 

Obviously,  each  such  solution  describes  a  certain  stationary,  horizontal 
vortidty  notion.  From  (12)  it  is  evident  that  each  stationary  solution  can  be 
determined  completely  by  assigning  one  function  v8(x,  y,  z). 

It  1s  not  hard  to  se#  that  any  solution  for  which  9  2  0  is  stationary. 
Actually,  if  9  s  0,  it  follows  from  the  first  two  equations  of  (9)  that  ♦  and 
p  arc  mot  functions  of  t  (and  are  connected  by  the  relationship  p  ■  <9).  Then 
the  third  equation  of  (9)  assumes  the  form 

c2  f*  ♦  px  -  0 

os 

end  the  solution  of  y  for  this  equation  with  boundary  conditions  (11)  is 
identical  with  zero.  From  the  last  two  equations  of  (9)  it  follows  that  p  is 
not  a  function  of  t  and  is  connected  with  p  by  the  relationship 


§?*r*p 

dz 


)  The  question  of  the  formula tlcn  of  the  condition  at  infinity  deserves- a  more 
detailed  examination.  Let  us  note  that  when  analyzing  wave  processes  it  turn*' 
out  that  the  condition  pw  — »0  is  too  weak  and  should  b#  replaced  fey  the  stronger 
condition  w^/p-vO  when  s— ioo. 


Tho  solution  of  equations  (9)  with  initial  data  (10)  will  be  stationary 
when,  and  only  when,  these  initial  date  satisfy  the  rel  atienship 

?.~o,  x.=o.  c**3— v^r-  <i^) 

It  follows  from  (12)  that  this  condition  is  necessary.  To  shear  that  it 
suffices,  lei.  us  differentiate  the  second  and  fourth  equations  of  system  (9) 
toth  respect  to  t  and  in  the  first  part  of  the  obtained  equations  e.colude  the 
derivatives  d(/dt,  dp/dt  and  dp/dt,  using  the  remaining  equations  of  (9).  As  a 
result  tro  get  the  following  systaas  of  fourth-order  time  equations  for  *  and  gi 

(-£  +  '*)’  +  +  e'** 

x  775-  “  V«  (fa  +  e’  *  +  £‘A?)  +  e  (  *' +  A9)  ’  <1S> 

Ths  solutions  of  9  and  y  for  these  equations  with  identical  boundary  con¬ 
ditions  (11)  is  determined  ccapletcly  by  the  following  initial  data,  '.then 
l  «■  0:  9  —  9»;  *■  “  ?*•  X  “X»-  *  "■  ”  ("ST  **“*)'  (16) 

Mien  conditions  (Hi)  are  fulfilled,  these  initial  data  revert  to  pero  end 
■consequently,  equations  (15)  have  only  the  trivial  solution  9  a  0  and  jsO, 
but  from  the  condition  9  5  0  it  follows  that  the  solution  of  equations  (9)  in 
this  case  trill  be  stationary. 

If  tho  initial  data  (10)  do  not  satisfy  relationships  (Hi),  the  solution 
of  equations  (9)  with  these  initial  dsti  can  be  represented  in  the  form  of  the 
sub  of  the  stationary  solution  which  can  be  determined  by  some  function  Vs 
and  of  the  non-statlonaiy  solution  of  9*,  0?,  p',  x'i  and  p'  which  satisfies 
the  initial  data.  When  ' 

t.O:  9.;  />'-*  —  x’-x»: 

In  order  to  find  the  function  jr3  which  determines  the  stationary  solution 
in  this  case,  we  nay  use  the  fact  that  equations  (9)  allow  an  invariant,  i.e., 
sone  expression  which  does  not  chango  with  time  and  which  is  a  linear  function 
of  the  unknown  functions  9,  p,  y,  and  p  .  let  us  use  the  folio-ring  method 
to  find  this  invariant.  I et  us  introduce  temporarily  the  new  notations  for  the 
unknown  functions 


'  Let  us  note  that  equations  (9)  also  permit  a  cer*ain  integral  invariant 
which  Is  quadratic  with  respect  to  the  unknown  functions  (the  energy  integral), 
/ctually,  using  equations  (5)  we  can  convince  ourselves  that  tho  following 
equality  holds«  .  ......... 

+ £ [*+**-  Ml  -  -(* ■ ♦*♦*)• 

from  which  it  follows  that  the  magnitude  [footnote  continued  next  page) 


yt  -  4>.  v«  -  ? .  v»  -  ?>  v»  -  x*  v»  ~  > 

end  writ*  equation*  (9)  in  the  form 

%-iU ,y.  (<-1.2 . 5). 

W  i-i 

where  Ajjj  are,  generally  speaking,  the  differential  operators 
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7e  will  seek  the  invariant  in  the  fern 


4-  intVi. 


i-i 


(21) 


where  yj»  generally  speaking,  are  linear  operators  independent  of  tin*  t.  Since 
djj/dt  should  eaual 

i-i  i-l  \t-l  / 

the  operators  n  satisfy  the  system  of  equations 


connected  with  the  system 


2l,ilttes0, 


^jAtiVt  —  0. 


(22) 


which  determines  the  stationary  solution  y^  -  (y^)8. 

It  is  easy  to  solve  equations  (22)  with  respect  to  the  operators  in  the 
given  case.  From  the  structure  of  the  first  and  fifth  columns  of  matrix  (20) 
it  is  evident  that  we  should  set  y2  •  ■  0.  In  this  case,  the  equation  ob¬ 

tained  using  the  third  column  is  identically  satisfied.  Thus,  only  the 


[/'+  dr 

is  independent  of  tame,  if  the  total  work  of  the  pressure  forces  on  the  sides 
of  a  cylindrical  container  of  base  S  is  equal  to  zero.  The  expression  in  the 
braces  is  the  density  of  the  complete  energy  of  the  disturbances,  and  the 
second  component  in  this  expression  is  the  density  of  potential  energy.  Thus, 
if  the  complete  energy  of  the  entire  field  of  disturbances  is  finite.  It  does 
not  change  with  time. 


1 


-  xo  - 


operators  n»  ****  My  differ  from  zero.  For  those  operators,  using  the 
second  and  fourth  columns  of  matrix  (20),  we  get  the  equations 

til  +  T.«'d  +  T»a  “  °.  T*  (?  +  *'  If)  +  T»  jV  “ 

Setting  yi  ■  A  wd  seeking  y3  and  Yj  in  the  form  of  linear  functions  of 
the  operator  d/dt,  tie  get 

ta-'if.  • 

Substituting  these  values  of  the  operators  Yj_  in  (21)  end  returning  to 
the  old  variables,  -we  get  the  follovrlng  invariant! 

(23) 

This  invariant  can  be  called  the  potential  vorticity  [1,  1*].  2,et  us  note 
that  from  the  third  and  fifth  equations  of  (9)  we  get  the  equality 

j 

"■'rlting  this  equality  with  z  -  0  and  remembering,  according  to  (11),  that 
the  magnitude  x  reverts  to  zero  when  z  •  0,  we  get  the  additional  invariant 

y.-y-cV.  (2U) 

where  the  asterisks  indicate  values  of  the  functions  when  z  »  0.  Let  us  note 
that  and  J2  do  not  contain  the  parameter  X,  i.e.,  they  do  not  change  fora 
when  transferring  to  e  quasistatic  approximation. 

Let  us  examine  the  case  where  the  initial  conditions  (10)  satisfy  con¬ 
ditions  (lL)  eve ryr.here  with  the  exception  of  a  certain  limited  region  of 
space  T.  The  non-stationary  component  of  the  solution  of  equations  (9)  under 
these  initial  conditions  will  be  of  the  nature  of  waves,  propagating  with  a 
certain  finite  velocity  from  region  V,  and  the  entire  energy  of  these  waves 
will  be  finite.  Since  this  energy  will  be  propagated  throughout  an  ever  in¬ 
creasing  volume,  the  values  of  functions  </*,  o',  p',  x'  and  p',  which  describe 
the  non-stationary  component,  will  tend  toward  zero  when  t  approaches  infinity, 
at  every  fixed  point  in  space,  and  in  this  case  a  complete  solution  of  equa¬ 
tions  (9)  will  tend  toward  their  stationary  component.  From  this  it  follows 
that  the  invariants  and  Jg  will  coincide  with  their  values  for  the  station¬ 
ary  components  of  the  solution  and  accordingly,  for  the  non-stationary  com¬ 
ponent,  those  invariants  are  equal  to  zero.  The  solutions  of  equations  (9)  for 
which  5  0  will  be  called  wave  solutions. 

Thus,  the  solution  of  equations  (9)  with  arbitrary  initial  data  (10)  can 
be  represented  in  the  form  of  e  sum  a)  of  the  stationary  solution,  for  which 
the  Invariants  and  J?  are  determined  according  to  the  initial  data,  and  b) 
of  the  wave  solution  which  satisfies  the  Initial  conditions  (17).  It  remains 
only  to  show  hem  we  find  the  function  which  determines  the  stetionaxy 


solution,  knowing  tha  values  of  the  invariants  and  Jj.  For  this,  let  us 


Here  we  get 

*’  _  j 

(25) 

(26) 

s(y  +  v  + 


Thus,  the  function  can  be  found  as  the  solution  of  an  inhomogeneous 
elliptical  equation  (25)  with  the  inhomogeneous  boundary  condition  (26)  (as 
the  second  boundary  condition  it  is  sufficient  to  require  that  be  limited 
as  s  approaches  Infinity).  In  order  to  avoid  inhomogeneity  in  the  boundary 
condition  here,  we  may  write  the  invariants  and  Jj  in  a  single  integral 
fora.  Thus,  let  us  examine  the  following  composite  invariant 

-  j  y.rf* + + £  5  8f.fr,*  - 

•  ••  •  f  •  *  ■  . 

-  ....  -  (27) 

-  J  M  -  l?)  *  -  •'  J  '-y2  -  $  (A*  -  rfr. 

This  invariant  is  the  integral  form  of  the  potential  vorticity. 

In  the  quasistatic  approximation,  by  excluding  the  magnitude  p-  -  i  2S 

e  ds  . 

from  expression  (27)  and  by  switching  from  the  integration  variable  s  to  the' 
new  variable  X  -  F/P*,  we  can  give  (27)  in  the  fomTK^TKOT  EEPKOBU' 

*)-'§•  ;•  <a) 

where  • 


is  the  non-dimensional  stability  parameter.  A.  S.  Monin  [It]  found  the  in¬ 
variant  (26)  for  a  quasistatic  barocllnic  atmosphere. 

Expressing  invariant  (27)  by  using  formulas  (12),  and  setting  V'j/p 
•XV,  we  get  the  folloring  equation  with  respect  to'pi 


(28) 

«  r  *- 1  h,-y\  t 
~  i  i *  »  V  t.  /  r* 

(29) 

*  Replacing  in  this  equation  the  magnitude  e  /<J  by  its  mean  value  in  the 
troposphere,  the  components  with  the  previous  derivatives  can  be  written  in 

,2  2 

the  form  Aya  ♦  — ^  ,  from  which  it  is  evident  that  the  characteristic 

scales  of  the  synoptic  processes  along  ths  vertical  and  ths  horitonsal  art 
most  naturally  determined  such  that  their  ratio  is  ot//fi  ,  which  is  approxi¬ 
mately  one  one-hundredth. 
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where  tho  right  side  can  be  determined  from  the  initial  data,  using  formula 
(27),  The  solution  of  this  equation  (in  the  caso  a2  is  const)  is  given  in  [I*], 

i 

» 

3.  THE  V.'AVE  SOLUTION 

.liter  finding  the  function  ys  which  determines  the  stationary  solution, 
from  initial  conditions  (10)  we  can  separate  the  part  [the  initial  conditions 
(17)]  which  determines  the  nave  solution,  '."hen  finding  this  solution,  one  of 
equations  (9)  can  be  replaced  by  the  condition  that  the  potential  vertieity 
(27)  be  equal  to  zero,  so  that  the  aquations  which  describe  the  wave  solution 
are  fourth-order  time  equations.  In  principle,  these  equations  make  it  pos¬ 
sible  to  express  f ,  p',  and  p*  by  0  and  x»  without  integrating  with  respeot 
to  time,  while  the  magnitudes  9  and  x  are  determined  from  equations  (15)  with 
initial  conditions  (16),  not  containing  function  Va.  Therefore,  it  is  not 
necessary  to  know  the  stationary  solution  to  find  the  wave  solution. 

Sin ce  the  coefficients  of  equations  (15)  are  independent  of  x,  y  end  t, 
these  equations  have  partial  solutions  in  the  form  of  harmonic  waves  with  , 

amplitudes  which  are  functions  of  z 

V  (*,  y,  t,  1)  -  O  (i)  «*  <*••+*•*  •*', 

X(r,  t)  =»  X  (31) 

where  k^  and  kj  are  the  horizontol  wave  nuabers  which  can  be  arbitrary,  and  0 
the  frequencies  to  be  determined.  Unce  the  dynamic  equations  which  we  are 
examining  are  linear,  and  we  are  not  considering  either  the  influx  or  dissi¬ 
pation  of  energy  in  them,  the  energy  of  waves  of  type  (31)  cannot  change  with 
time.  Accordingly,  frequencies  a  must  be  real. 

The  solution  of  equations  (15)  with  arbitrary  initial  data  can  be  repre¬ 
sented  in  the  form  of  the  suoorposition  of  elementary  wave  solutions  (31)  with 
all  possible  values  of  kj_  and  kg.  Therefore,  the  explanation  of  the  nature  of 
tho  elementary  naves  (31)  makes  it  possible  to  draw  definito  conclusions  as  to 
the  nature  of  any  nave  solution  of  the  exanined  equations  of  hydrodynamics. 

In  order  to  otudy  waves  (31)  we  have  only  to  find  their  amplitudes  S(z)  and 
X(z)  and  determine  the  frequencies  o  (i.e.,  find  the  frequency  spectrum  of  the 
elementary  wave  aolutlona).  ‘Substituting  functions  (31)  in  equations  (15)  we 


o 
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get  the  following  equations  for  amplitudes  S(z)  and  X(z)i 

(i>  +  *•«•  — #*)  <t>  f*-V  -H  «*  -37" 1 
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Let  us  note  first  that  these  equations  have  a  non-trivial  partial  solution 
in  which  X  S  0  (since  this  solution  satisfies  tho  boundary  conditions  it  has 
a  physical  significance).  The  first  equation  of  (32)  sheers  that  the  frequen¬ 
ce*  corresponding  to  this  solution  are  determined  by  the  formula 

,«-*•  +  **<*.  (33) 

According  to  the  second  equation  of  (32)  the  amplitude  B(s)  should  satisfy 
the  aquation  ‘ 


so  that  it  has  the  fora 


-sT$ 

®  (t)  -®*«  •  -  ©•  (p/py*. 


(31*) 


Evidently,  the  waves  corresponding  to  this  solution  embrace  the  entire 
atmosphere  at  onoe,  propagate  only  horizontally  and  do  not- cause  vertical 
variations  of  air  particles.  In  the  quasistatic  approximation  (X  »  0),  these 
wares  retain  their  fora.  The  maxi  nun  group  velocity  of  motion  of  packets  of 
similar  waves  (the  velocity  of  the  front  of  the  wave)  is  equal  to  the  speed  of 
sound  c.  The  indicated  waves  are  completely  analogous  to  the  waves  in  a  quasi¬ 
static  barotropic  atmosphere  studied  by  A.  Obukhov  (1)  when  he  examined  the 
question  of  the  adaptation  of  hydrodynamic  fields.  Such  waves  were  singled 
out  by  A.  S.  Honin  [!*]  for  the  esse  of  a  quasistatic  baroclinic  atmosphere. 

In  the  future,  we  will  call  these  waves  two-dimensional. 

Turning  to  a  study  of  the  solutions  for  which  X  «£  0,  for  purposes  of  sim¬ 
plifying  the  computations  we  will  limit  ourselves  to  an  examination  of  the  case 
where  the  coefficients  c2  and  ?  in  equations  (32)  can  be  considered  constant 
[this  condition  is  accurate  for  an  isothermal  atmosphere  in  which  0^  ■  xgH  and 
p  ■  (x-l)g,  where  H  a  p*/gp*  is  the  height  of  the  unifora  atmosphere].  In  this 
case,  each  of  the  functions  9  and  X  satisfies  the  following  equation,  which 
stems  from  (32) 

<-->(£  +  05) 

This  equation  has  a  partial  solution  in  the  fora  [exp  (-I>ia)z],  whqre. 

V  and  m  ara  real  numbers.  .  Substituting  this  function  in  (35)*  we  get 
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ll*-.*)  [(-  M  +  /«)*  +  "  +  <«*>  +  Vr]  -  ($  “  x,%)  “  °* . 

Reseabering  that  a  is  real  and  equating  the  imaginary  part  of  the  obtained 
equality  to  taro,  we  get  U  •  (p*g)/2o2.  Using  this  result  and  equating  the 
real  part  of  the  equality  to  zero,  we  get  the  following  relationship! 

+  -vj  _  it(t 


Let  us  note  that  functions  ,  and  x  of  the  fora 


f.  X~*  *  • 


correspond  to  the  examined  partial  solution  of  equation  (35),  plane  has** 

nonic  naves  with  horlzonUi  wave  nrabers  ^  and  k2  and  a  vertical  nave  nuaber  . 
b,  nhcse  amplitudes  decrease  exponentially  with  height.  These  wavbj  have 
physical  significance  only  when  m  /  0,  since  when  o  ■  0,  no  combination  of 
such  naves  satisfies  the  boundary  condition  x — »0  (z— -jC).  Accordingly,  the 
propagation  rate  of  such  waves  always  has  a  vertical  component,  i.e.,  these 
waves  are  essentially  three-dimensional.  He  will  call  such  w eyes  ^sternal 


waves. 

Internal  waves  with  any  wave  numbers  can  enter  into  the  solution  of  the 
Cauchy  problem  for  equations  (15)  with  arbitrary  initial  data.  Therefore,  in 
(36)  the  magnitudes  k  and  a  must  be  considered  arbitrary  and  then  (36)  is  an 
equrtion  with  respect  to  o,  noting  that  the  left-hand  part  of  equality  (36) 
is  not  negative  and  considering  that  •  „ 

I- <*<*#7 

we  can  easily  convince  ourselves  that  the  roots  of  o2  in  equation  (36)  can  be 
found  only  within  the  intervale 

*•<••<•£.  (38) 

whereupon  (due  to  the  fact  that  the  roots  of  o2  aro  continuous  functions  of  lc? 
end  b2),  all  points  of  these  intervals  are  possible  frequencies  of  internal 


waves.  Thus,  two  fardlies  of  internal  waves  can  occur  in  the  atmosphero,  i.e., 
waves  with  frequencies  freo  the  first  or  from  the  second  interval  of  (38).  To 
analyze  the  nature  of  these  waves,  let  us  write  the  solutions,  of  equation  (36) 
for  the  square  of  the  freqrency  o?  in  the  form 
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where  L  •  e/4  is  the  character!  s  tic  scale  of  the  horizontal  movements  of  a 
compressible  medium  in  a  Coriolis  force  field,  first  introduced  in  [1], 

From  (39)  it  is  evident  that  the  influence  of  the  Coriolis  force  is 
essential  only  with  very  small  horizontal  wave  numbers  k  (l.e,,  for  waves 
whose  horizontal  scales  are  not  small  compared  with  L).  In  order  to  explain 
the  nature  of  waves  with  frequencies  ofl  and  let  us  first  digress  from  the 
effect  of  the  Coriolis  force  (i.e.,  let  us  sot  l  -  0  and  accordingly,  L  -  oo). 

9 

Let  us  assume  that  X  •  1  and  examine  the  isothermal  atmosphere  in  which  o  ■ 

•  xgH  and  0  •  (x-l)g,  Here  formulae  (39)  assume  the  fom 
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Let  us  remember  that  we  are  examining  polytropic  processes  during  which 
the  magnitude  pp-^A  is  retained  in  liquid  particles..  In  particular,  when 
x  -  oo,  the  value  retained  will  be  density  p,  so  that  this  case  corresponds  to 
a  non-comp resslble  modium.  V/hen  x  — i  oo,  all  frequencies  revert  to  infinity, 
i.e.,  the  waves  with  these  frequencies  disappear.  But  .in  a  non-compreeeible 
medium,  only  variations  with  frequencies  o^  may  arise,  determined  from  (ItO) 
within  the  boundazy  x  — >  oo  : 

-t  |  H 

*'  +  «*  +  ijjl 

Evidently,  waves  with  these  frequencies  are  internal  gravitational  waves. 
A  description  of  them  can  be  found,  e.g.,  in  Prandtl's  book  [8].  Further, 
when  x  ■  1,  the  value  which  is  maintained  will  be  temperature,  so  that  this 
.  case  corresponds  to  isothermal  processes,  with  respect  to  which  isothermal 
stratification  is  neutral.  V/hen  x— }1,  all  frequencies  Og  revert  to  zero,  so 
that  variations  with  these  frequencies  disappear  (they  convert  in^o  stationary 
yorticity  motions).  Thus,  with  isothermal  processes  in  an  isothermal  medium, 
only  variations  with  frequencies  o4  may  occur,  which  can  be  determined  from 
(1*0)  within  the  limit  x — ,1« 


-  iO  • 
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Evidently,  waves  with  these  frequencies  are  acoustic.  An  analogous  dis¬ 
cussion  shows  that  two-dimensional  waves  with  frequencies  (33)  are  acoustic. 

Returning  to  the  case  of  adiabatic  processes  during  which  x  ssl.li,  we  say 
conclude  that  the  frequencies  °a  dote  rained  by  formulas  (liO)  correspond  to 
ecoustic  naves,  distorted  due  to  the  absence  of  isothursy  in  neutral  strat¬ 
ification,  while  frequencies  og  correspond  to  Gravitational  naves  distorted 
because  of  the  compressibility  of  the  air.  '  An  analogous  interpretation  can  be 
applied  to  naves  corresponding  to  frequencies  (39)  end  containing  an  additional 
distortion  due  to  the  effect  of  the  Coriolis  force. 

Assuming  X  •  1  in  formulas  (39)  *nd  making  use  of  the  fact  that  with  suf¬ 
ficiently  large  k  and  m,  the  subtrahend  under  the  radical  in  these  fonaules  is 
small  (so  that  wo  may  assume  approximately  that  ^  1  -  J.  ~1  -  i),  ws  get  the 
following  approximate  formulas  (which  are  the  moro  accurate,  the  larger  ere  k 
and  m)t 
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2  2 

The  inequalities  shown  here  for  oa  and  og  are  strict  with  any  k  and  m. 
Thus,  naves  nith  frequencies  from  the  first  interval  of  (38)  are  gravitational, 
those  from  the  second  interval  are  acoustic.  In  the  case  of  an  isothermal 
atmosphere,  setting  k  -  1.1  and  H  ■  Q  km,  we  get  the  following  numerical  ine¬ 
qualities  for  the  periods  of  acoustic  end  gravitational  waves: 


r.- 


*"<300  sec;  Tt  -  •*"  >  330  see 

*•  i 


In  an  atmosphere  in  which  temperature  drops  with  height,  the  interval 
between  Ta  and  Tg  increases. 

The  frequencies  of  slight  variations,  possible  when  using  the  quaslstatle 
approximation,  are  obtained  from  frequencies  (39)  within  the  limit  X— *0. 

Using  this  limited  transition,  all  frequencies  og  revert  to  infinity,  i.e., 
the  internal  acoustic  waves  are  completely  filtered  out.  Here  the  group  ve¬ 
locity  of  the  vertical  propagation  of  packets  of  internal  acoustic  waves 
doa/dm  also  reverts  to  infinity,  so  that  in  the  quasistatic  approximation  the 
propagation  rate  of  disturbances  along  the  vertical  appears  infinitely  large. 
Uhan  the  quaslstatle  approximation  is  used,  the  only  Internal  waves  possible 
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are  the  gravitational  on* a  whos*  frequencies  are  determined  by  the  fomula 


Comparing  this  formula  with  the  corresponding  formula  (J*3),  we  see  that 
the  frequencies  of  gravitational  waves  in  the  quasistatio  approximation  are 
too  high,  snd  they  are  the  greater,  the  larger  the  horizontal  wave  number  k;  ' 
however,  for  the  -case  k  ^  a  (long  waves),  this  change  of  a  remains  insignifi¬ 
cant. 

Visaing  up  the  results  of  this  analysis,  we  see  that  the  slight  variations 
which  can  occur  in  the  atmosphere  with  arbitrary  initial  data  break  down  into 
three  classes*  two-dimensional  waves,  internal  acoustic  waves  and  internal 
gravitational  waves. 

The  role  of  the  quasistatio  approximation  basically  reduces  to  the  fil¬ 
tering  out  of  internal  acoustic  waves  (a  secondary  effect  is  a  certain  dis¬ 
tortion  of  internal  gravitational  waves).  Therefore,  when  the  internal  acous¬ 
tic  waves,  which  arise  due  to  the  disruption  of  static  equilibrium  A  the 
initial  moment  of  time  in  a  certain  region  of  space  V,  scatter,  we  may  consider 
that  the  atmosphere  has  basically  already  adjusted  itself  to  the  state  of 
static  equilibrium.  The  time  of  adaptation  to  the  state  of  static  equilibrium 
•Is  of  the  same  order  of  magnitude  as  the  time  it  takes  for  the  front  of  the 
internal  acoustic  waves  to  traverse  the  main  part  of  the  atmosphere.  Since 
this  front  mores  at  the  speed  of  sound  c,  the  Indicated  time  is  only  several 
minutes  (in  one  minute  a  sound  wave  traverses  a  layer  20  km  thick).  After 
static  equilibrium  has  been  attained,  the  process  of  the  adaptation  of  the 
atmosphere  to  the  state  of  geostrophic  equilibrium  continues,  whereupon,  on  an 
average,  such  a  state  is  established  throughout  the  atmosphere  after  the 
scattering  of  two -dimensional  waves,  while  somewhat  later,  after  the  scattering 
of  slow  (see  [li])  internal  gravitational  waves,  geostrophic  equilibrium  Is  es¬ 
tablished  at  all  altitudes.  The  rate  of  this  latter  process  depends  essen¬ 
tially  on  atmospherio  stratification. 


Academy  of  Sciences  of  the  USSR 
Institute  of  Atmospherio  Physio* 


Submitted 
.  10  April  1953 


r 


0 


O 


-  18  - 


LITERATURE  CITED 


[1]  Obukliov,  A.  IS.  "K  voprosu  o  ceostroficheskom  vetre"  (The  question  of  the 

Geostrophic  wind),  Akadeniia  Kauk  SSSS,  Leningrad,  Izvestilg.  Serila 
Geottrafichoskala  1  Geofizicheskaia.  13(4) «  261-306.  1949. 

[2]  Charney,  J,  G.  "On  the  scale  of  ataospheric  notions,"  Goofysi'ske 

Fubiika.rloner.  17(2),  1948. 

[3]  Kibel',  I.  "0  prisposoblenii  dvizi.eniia  vozdukha  k  geostroficheskouu* 

(The  adactation  of  air  movement  to  the  gecstrophic),  Akadeniia  Kauk 
SUSR,  Leningrad,  Dokladv,  104(1):  60-63,  1955. 

[4]  Konin,  A.  S.  "Izmcneniia  davleniia  v  baroklinnoi  ataosfere"  (Pressure 

changes  in  a  baroclinic  atmosphere),  Akadeniia  Kauk  S3SR,  Leningrad, 
Izvertlia,  Serila  Geoflsicheskala,  IioT"E»  4?7-5l4, 

[5]  Dorodnitsyn,  A.  A.  "Vozaushchenie  vozdushnogo  potoka,  vyzvrnnoe  nerov- 

nost'iu  noverkhnosti  ZerJi"  (Tne  disturbance  of  an  air  streaa  caused 
by  the  unevenness  of  the  earth's  surface),  Olavnaia  Geofizlcheskaia 
Observatorila,  Trudy,  23(6),  1938. 

[6]  Lyra,  F.  "fiber  den  Einfluss  von  3odenerhebuncen  auf  die  Strflnung  einer 

stabil  geschichteten  Atuosphslre"  (The  effect  of  ground  elevations  on', 
the  current  of  a  stable  stratified  atrosphere),  -3eitr3ge  zur  Physik 
der  frelen  Atnosoharo.  26(4):  193-206,  1940. 

[7]  Iaglcn,  A.  K.  "Dinsnika  krupnoisanshtabnykh  protscssov  v  barotropnoi 

atnosfere"  (The  dynamics  of  la  ryj- scale  processes  in  a  barotropic 
atmosphere ),  Akadeniia  Kauk  S  '..enin-rsd,  Izvestii?,  Serila  Geo¬ 
fizlcheskaia,  Ho.  346-369.  1953^ 

[8]  Prandtl,  L.  Gldroaeronekhanika  (Hydroaeronechanics).  ilosccw,  Izdatell- 

stvo  Znostramoi  Xiteratury,  1951. 


